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Abstract 



We consider effects of the isocurvature perturbation in the framework of the 
tracker-type quintessence models. During the inflation, fluctuations in the ampli- 
tude of the quintessence field are generated, which provide isocurvature component 
of the cosmic density perturbation. Contrary to the conventional notion, it is shown 
that effects of the isocurvature fluctuation may become sizable in some cases, and in 
particular, the cosmic microwave background angular power spectrum may be signif- 
icantly enhanced due to the effect of the isocurvature mode. Such an enhancement 
may be detectable in on-going and future experiments. 



Recent cosmological observations suggest that there exists a dark energy which must 
be added to the matter density in order to reach the critical density [1]. Although the 
cosmological constant is usually assumed as the dark energy, in the past years, a slowly 
evolving scalar field, dubbed as "quintessence" has been proposed as the dark energy 0. 
There are some differences between the cosmological constant and quintessence. Firstly, 
for the quintessence, the equation-of-state parameter ojq = Pq/pq varies with time, whilst 
for the cosmological constant, it remains a fixed value uj^ = —1. Secondly, since the 
quintessence is a scalar field, it fluctuates. 

One of the observational discriminator of these models is the cosmic microwave back- 
ground (CMB) anisotropies. Many authors have studied the effects of the quintessence field 
on the CMB anisotropy in various models. As mentioned above, the quintessence field fluc- 
tuates, so its fluctuations should be taken into account in calculating the CMB anisotropy. 
Generally, the initial fluctuations of the quintessence fields are generated during inflation. 
Since these fluctuations behave as the isocurvature mode, there can exist both adiabatic 
and isocurvature perturbations. Implications of the isocurvature modes have been studied 
in [HI, and in particular for the cosine-type quintessence model, it was pointed out that 
the isocurvature fluctuations may be detectable in on-going and future experiments by the 
present authors 

In this letter, we focus on the effects of the isocurvature fluctuations on the CMB 
anisotropy in the tracker-type quintessence models. The tracker-type quintessence models 
have attractor-like solutions which alleviate the initial condition problems. When we take 
initial conditions such that the quintessence field starts to roll down the potential at early 
epoch, damping effect on fluctuation of the quintessence field is significant. So, in this case, 
even if we take a nonzero value of primordial quintessence-field fluctuation, the fluctuation 
goes to zero because of the damping effects. But we will show that, in some cases, the 
damping effect is not so significant. In such cases, the effect of the isocurvature fluctuations 
may be sizable and affect the CMB angular power spectrum. 

Since the effects of the isocurvature fluctuations depend on models of quintessence, 
first, we classify the tracker-type quintessence models. Although many models have been 
proposed P,|3L H-EDl, they can be classified into two groups, using the evolution of the en- 
ergy density of the quintessence pq. The evolution of the energy density of the tracker-type 
models can be written in a simple form. Tracker-type quintessence models, as previously 
mentioned, have an attractor-like solution in a sense that a very wide range of initial condi- 
tions converge to a common evolutionary track. We call the epoch when this attractor-like 
solution is realized as "tracking regime." In the tracking regime, ujq is almost constant and 



1 — ujq is significantly different from zero |21[ . This means that the kinetic and potential 



energy of the quintessence have a fixed ratio. Thus, the relation Q 2 oc V(Q) holds in 
the tracking regime. When this relation is satisfied, we can show that the energy density 
of the quintessence is proportional to a~ n , where n is a constant. On the contrary, the 
energy density of the dominant component of the universe is written as po oc aT m (for the 
matter and radiation dominations, m = 3 and 4, respectively). When n = m, the energy 
density of the quintessence can be a significant fraction of the total energy density of the 
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universe from earlier epoch. On the contrary, when n < m, pq cannot be significant at 
early epoch.j^] For a purpose of this letter, we will consider two models as representative 
models; "AS model" || and "Ratra-Peebles model" [13] for models with n = m and n < m, 
respectively 

Before we show the evolution of the quintessence energy density from numerical calcu- 
lations, we investigate it analytically. To study the evolution of the quintessence fields, we 
write down the basic equations. In a spatially flat homogeneous universe which contains a 
perfect fluid with energy p and pressure p and a quintessence field Q with energy pq and 
pressure pq, the Friedmann equation becomes 

# 2 = ^(p + pq), (i) 

where H is the Hubble parameter and M* ~ 2.4 x 10 18 GeV is the reduced Planck scale. 
The energy density and pressure of the quintessence field are 

Pq = \q 2 + V(Q), P q=Iq*-V(Q), (2) 

where the dot represents the derivative with respect to time t. The quintessence field Q 
obeys the equation of motion 

Q + 3HQ + % = 0. (3) 
dQ 

The energy conservation equation for each component is, in case where there is no energy 
exchange 

p = -3H(p + p). (4) 



Thus, when the universe is dominated by a perfect fluid with the energy density pr> oc a 
the scale factor behaves as 



-m 

j 



a oc t 2/m . (5) 



In this case, from Eqs.(@) and 



First, let us discuss the evolution of pq for the AS model which has the potential of the 
form 

V (Q) = \(Q - bf + a] e- XQ = f(Q)e' XQ , (7) 



The case with n > m is also possible [f22| , but almost all model which have been proposed so far, can 
be categorized into the type-I or the type-II models. So we will not consider the case with n > m in this 
letter. 
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where a, b and A are model parameters. Since f(Q) changes more slowly than e~ x ® in the 
tracking regime, in the analytic investigations below, we approximately treat f(Q) as a 
constant. Substituting Eq.(^) into Eq.(0), we obtain 

Q + — + Q - f(Q)^- XQ = o. (8) 

mt 

The solution to this equation is Q = AlnXBt, where A and B are constants. In the 
tracking regime, Q 2 is proportional to V(Q), thus pg oc t~ 2 in this model. Using Eq.(|5|), it 
is shown that 

Pq oc a~ m . (9) 

Therefore, in the AS model, the energy density of the quintessence has the same a depen- 
dence as that of the dominant component of the universe. This means that the energy 
density of the quintessence can be a sizable fraction of the total energy density of the 
universe at early epoch. In fact, the above result does not depend on the detailed struc- 
ture of the function f(Q) as far as f(Q) is a slowly varying function on the tracking 
regime. Thus, the same discussion is applicable to other types of quintessence models with 
V(Q) = f(Q)e- x Q. 

Next, we consider the Ratra- Peebles model which has the potential of the form, 

V{Q) = (10) 

where A and a(> 0) are model parameters. With this potential, the equation of motion of 
the quintessence field becomes 

6 • A 4+a 

This equation has the solution Q = Ct u where v = 2/(2 + a) and C is a constant. Since a 
is taken to be positive, the value of v is < v < 1. The energy density of the quintessence 
in this model becomes 

p Q oc a~ m{ ±-»\ (12) 

Thus, in the Ratra-Peebles model, pq decreases more slowly than that of the dominant 
component of the universe, which implies that the energy density of the quintessence cannot 
be a significant fraction of the energy density in the early universe. 

Now we show evolutions of the energy density of the quintessence models from the 
numerical calculations. In Fig.|l|, we show the evolution of the quintessence energy density 
for the AS model, where we take several values of the initial amplitude for the quintessence 
field. From the figure, we can read off when the quintessence field starts to enter the 
tracking regime for several values of the initial amplitudes. We can also see the property of 
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the AS model; the energy density of the quintessence field can have O(0. 1) contributions 
to the total energy density of the universe from early epoch. 

In Fig. [| the evolution of the energy density are shown for the Ratra-Peebles model. 
Different from the AS model, the quintessence energy density of this type of model is 
negligible for z 3> 0(1). This fact has an important implication to the effects of the 
isocurvature fluctuations on the CMB angular power spectrum. 

Next let us discuss the effects of the isocurvature fluctuations on the CMB angular 
power spectrum. For this purpose, we study the evolutions of the fluctuations of the 
quintessence analytically. We decompose the quintessence field Q as 

Q(t,x) =Q(t) +q(t,x), (13) 

where q is the perturbation of the amplitude of the quintessence field. The equation of 
motion for q is, in the conformal Newtonian gauge, 

fa\- 2 n 9 d 2 V(Q) ,i - x =- f a\ 2 dV , , 

where the perturbed line element in the conformal Newtonian gauge is given by 

^ 2 

Mo 

—) + 2^)dr 2 + (1 + 2<f>)5 ij dx i dx j ] , (15) 

a Q J 1 J 

where a is the scale factor at time t, a® the scale factor at the present time and r is the 
conformal time. 

From Eq.fll|]), we can show the damping behavior of the quintessence field in the 
tracking regime. To solve this equation, we write down the second derivative of the potential 
as a function of the sound speed c 2 of a quintessence field Q H 



ds 2 = -(1 + 2V)dt 2 + ( — ) (1 + 2§)5 ij dx l dx 3 



H 3 , 9 „ « 



dQ 2 2 s 2 
where the sound speed of the quintessence field is written as 



(16) 



2 PQ Q - dV/ dQ 



_^ = ? (17) 
Pq Q + dV/dQ 

Since the kinetic energy is proportional to the potential energy in the tracking regime, the 
relation Q oc dV/dQ holds. So the sound speed of the quintessence field c 2 is a constant 
during this regime. Therefore c 2 can be set to zero in Eq. fllCf) . When the quintessence is 
a subdominant component of the universe, the Hubble parameter H can be written using 
the equation-of-state parameter ljjj = Pd/pd of the dominant component of the universe, 

H=— — -. (18) 

3(l + u D )t v ; 
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To study the isocurvature mode at the superhorizon scale, we neglect the k dependence and 
the right hand side of Eq.([l4]). Then, using the above relations, the equation of motions 
can be written as 

« + n 1 , u g + n X ~ ^ f2 ( 1 - c ')( c ' + 2 )9 = 0- ( 19 ) 

This equation has power-law solutions like q oc where the power-law index £ is 
wo - 1 



2(u^ + V 



(20) 



Since is equal to ujq in the tracking regime and ojd,q cannot be larger than unity, 
the real part of £ is always negative. Therefore q damps with time. It follows that if we 
take q ^ initially, q damps to zero in case that the quintessence field experiences long 
period of tracking But, in the case that the quintessence field enters the tracking 

regime at later time, the fluctuations may not damp so much. In this case, the isocurvature 
fluctuations (i.e., q ^ 0) can affect CMB anisotropics. 

The primordial fluctuations in the quintessence amplitude is generated in the early 
universe, probably during the inflation. When the mass of the quintessence field is negligible 
compared to Hubble parameter during inflation H in f, the primordial fluctuation is 



= ff ■ (21) 

To parameterize the size of the isocurvature contribution, we define the following quan- 
tity r q : 

- = _J*5_ (22) 



where q^ is the fluctuation of the quintessence field in the deep radiation-dominated epoch 
when the quintessence field is slow- rolling, while ^jj^ is \I> from the adiabatic mode (i.e., 
contribution from the inflaton-field fluctuation) in the deep radiation-dominated epoch. 
For a given model of slow-roll inflation, is 



T (adi) _ 4 ( ff inf Vjtf , 

^ " 9 v^rw^) ' (23) 

where is the inflaton potential and V{ ni is its derivative with respect to the inflaton 
field x- 



# 2 If m q ^ H[ n { , the primordial fluctuation damps to zero rapidly and cannot affect the CMB anisotropics. 



So we do not consider such a case 
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The parameter r q depends on the detailed scenario of cosmology. In a simple case with 

<?RD = Qmii 

For example, for the chaotic inflation with Vi n f oc x p where p is an integer, r q is given by 

9pM, 



q | chaotic 



4 X (h 



COBEj 



(25) 



where x(^cobe) represents the inflaton amplitude at the time when the COBE scale crosses 
the horizon. Numerically, we found r q ~ 0.3 — 0.6 for p = 2 — 10. The ratio r q may 
become larger, however, in more complicated cases. For example, if the kinetic function 
for the quintessence field Zq varies during inflation, q also changes. This may happen if 
Zq depends on the inflaton field. Therefore, in general the ratio r q is model dependent 
and we treat free parameter in our analysis. 

Now we show the results from numerical calculations. We used the modified version of 



CMBFAST [25] to calculate the CMB angular power spectrum C; which is defined as 



(AT(x, 7)AT(f , f )) = A. ]T(2Z + 1)C z P z (t ■ f), (26) 

4vr , 

where AT(x, 7) is the temperature fluctuation of the CMB pointing to the direction 7, and 
Pi is the Legendre polynomial. The average is over the position x. 

In Fig. 0, we plot the CMB angular power spectrum C\ in the AS model with several 
values of the initial amplitude Q in and the ratio r q .f^\ For comparison, we also plotted C; 
in the cosmological constant case. 

First, we discuss the adiabatic case (i.e., r q = 0), in particular, focusing on differences 
between the quintessence and the cosmological constant cases. One can see that the lo- 
cations of the acoustic peaks may differ in two cases. In our case, locations of the peaks 
mostly depend on the angular diameter distance to the last scattering surface (LSS) r^r*), 



and the location of the n-th peak in the I space is estimated as |26|.|27 



l n ~ ^rnvr, (27) 

where r s (T*) is the sound horizon at the recombination. The angular diameter distance 
becomes smaller when the energy density of the universe is large since the expansion rate 
becomes larger in this case. The energy density of the universe in the earlier epoch becomes 
larger in the case with the quintessence. Thus, in the quintessence case, the locations of 
the acoustic peaks are shifted to lower multipoles due to this effect. Such a shift is observed 
in Fig. 

# 3 With the cosmological and model parameters used in Fig. ||, we checked that theoretical predictions 
of Ci are in good agreements with current observational data with proper normalization. 
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Another point we should mention is the height of acoustic peaks (relative to C\ at lower 
multipoles) . In the AS model, the driving effect is important to understand the height of 
the peaks; in the radiation-dominated epoch, the amplitudes of the acoustic oscillation is 
boosted just after the horizon crossing due to the decay of ^ . Compared to the cosmological 
constant case, the expansion rate of the universe becomes larger in the quintessence case, 
resulting in faster decay of the gravitational potential ^ '. This effect enhances the driving 
effect. Thus, in particular in the AS models, height of the acoustic peaks are enhanced 
relative to lower multipoles ||. 

Next we consider the case with r q ^ 0. For this purpose, it is instructive to study the 
behavior of the perturbations generated by the fluctuation of the quintessence amplitude. 
Evolutions of the perturbations are governed by the following equations; for the metric 
perturbation, 



fc 2 $ = 47rG - p to t 



a\ 2 r„ _ .V t 



tot 



5 to t + 3W(1 + 



(2f 



and for the photon, at the superhorizon scale, 



5; = -^V 7 - 4$', V? = *Qfi y + *) l (29) 

where 7i = a' /a, 8 1 and V 1 (<5 to t and Vtot) are density perturbation and velocity perturbation 
of photon (those of total matter density) which are defined in the Newtonian gauge, and 
the prime represents the derivative with respect to the conformal time r. In the deep 
radiation-dominated epoch when the quintessence field is slow-rolling, the fluctuations in 
energy density, pressure and velocity of the quintessence are approximated as 

dV ( a \ ~ 2 

5p Q ~ -5p Q ~ — q KB , (p Q + p Q )V Q = k i^—J Q'q RD . (30) 

In the deep radiation-dominated epoch when the scale we concern is out of horizon, we can 
study the behavior of the perturbations by expanding each variables with the conformal 
time r. Neglecting the shear stress from neutrino (i.e., using II tot = 0), ^ ~ — $. Thus we 
find 

$ ~ CW 4 , 5 7 ^ 4Cg RD r 4 , V 1 ~ ^Ckq RD r 5 , (31) 

with 

1 dV 1 

6 " "126 rfQMf rad0 ' (32) 

where p ra do is the energy density of radiation at the present epoch. Here, we keep the 
leading term in r. In the radiation-dominated epoch, r oc t l l 2 and all of these variables 
(i.e., <5 7 and Vy) grow with time for the isocurvature mode. In addition, we can see that, 
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at the superhorizon scale, \l/ has the same k dependence as ^rd- Thus, if q^D has no scale 
dependence, \I/ is scale-invariant. 

This behavior can be confirmed by numerical calculations. In Fig. [|, we show evolution 
of the gravitational potential generated by the fluctuation in the quintessence amplitude. 
When the quintessence is in the slow-roll regime, $ is proportional to a 4 . (Notice that, in 
the radiation-dominated epoch, a oc r.) At the time when the quintessence field enters the 
tracking regime, ^ ceases to grow since q starts to damp. Thus, ^ has its maximum value 
at this epoch. Consequently, the CMB angular power spectrum is the most enhanced at 
the angular scale which enters the horizon at the time when the tracking regime starts. 

In the AS model, since Qq can be O(0. 1) in the tracking regime, fluctuation of the 
quintessence field may significantly affect the angular power spectrum C\. To see the effect, 
in Fig. [], we plot the CMB angular power spectrum for the case with r q = 5. (Notice 
that the enhancement due to the isocurvature contribution, Ci{r q ^ 0)/C[(r q = 0) — 1, 
is proportional to r q .) As we can see from the figure, when r q ^ 0, the angular power 
spectrum Ci is enhanced. 

To study the enhancement in C\ by the isocurvature fluctuation, we plot the ratio 
Ci(r q = 5)/Ci(r q = 0) with several values of the initial amplitude in Fig. |. We can clearly 
see that the multipole with the largest enhancement depends on the initial amplitudes of 
the quintessence fields. As we discussed before, this angular scale is determined by the time 
when the tracking regime starts; the effect of the isocurvature fluctuation is not significant 
for the angular scale much larger than and much smaller than this scale because of the 
damping effect of q or smallness of Qq. As we can see, the later the quintessence field 
enters the tracking regime, C\ is most enhanced at lower multipole /. Since the epoch 
when the quintessence field enters the tracking regime depends on its initial amplitude, the 
enhancement in Ci may give information about the initial condition of the quintessence 
fields. In addition, with r q = 5, the enhancement in C\ can be as large as 50 % relative to 
the case with r q = when Q m = 121M*, which can be detectable in on-going and future 
satellite experiments Thus studying effects of isocurvature fluctuations provides 

us implications on model building of the quintessence, which is complementary to studying 
the adiabatic fluctuation. 

Finally, we comment on the case with the Ratra-Peebles model. We also calculated the 
CMB angular power spectrum in the Ratra-Peebles model. However, in this model, effects 
of isocurvature fluctuation are very small since Qq is very small in the early universe. For 
example, even if we take r q = 5 with the model parameters A = 5.122 x 10 6 GeV and 
a = 6, enhancement in C\ is less than a few %. 
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Figure 1: Evolution of the energy density in the AS model. The model parameters are 
taken to be A = 2.2/M^A = 0.01M, 2 and B = 123.6805M*. The initial amplitude are 
Q in = 123M* (solid line), 121M* (dotted line), and 118M* (dashed line). The cosmological 
parameters are taken to be f2 = 0.3, Vl b h 2 = 0.019 and h = 0.65, where h is the Hubble 
parameter in units of 100 km/s/Mpc. 
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Figure 2: Evolution of the energy density in the Ratra-Peebles model. The parameter we 
take here, A = 5.122 x 10 6 GeV and a = 6. The initial amplitudes of the quintessence 
field are Q in = (solid line), 0.5M* (dotted line) 0.1M* (dashed line). The cosmological 
parameters are the same as those in Fig.|l]. 
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Figure 3: The CMB angular power spectrum in the AS model. The initial amplitudes are 
(a) Q in = 123M*, (b) Q m = 121M*, and (c) Q in = 118M*. The cosmological and model 
parameters are the same as those in Fig. |1]. The overall normalization is arbitrary in this 
figure. 
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Figure 4: The evolutions of \^\ of purely isocurvature mode with (A) Qj n = 12 1M* and 
(B) Q in = 118M*. We take fc=lx l(T 4 /i Mpc" 1 (solid line) and k = 1 x l(T 2 /i Mpc -1 
(dotted line). The model and cosmological parameters are the same as those in Fig. |l|. 
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Figure 5: The ratio Ci(r q = 5)/Ci(r q = 0) in the AS model. The initial amplitude 
are taken to be 123M* (solid line), 121M* (dotted line), and 118M* (dashed line). The 
cosmological and model parameters are the same as those in Fig. [l|. Notice that the 
quantity Ci{r q ^ 0)/Ci(r q = 0) — 1 is proportional to r^. 
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